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Abstract. We study a one dimensional nonequilibrium lattice model with competing 
features of particle attraction and non-local hops. The system is similar to a zero 
range process (ZRP) with attractive particles but the particles can make both local 
and non-local hops. The length of the non-local hop is dependent on the occupancy 
of the chosen site and its probability is given by the parameter p. Our numerical 
results show that the system undergoes a phase transition from a condensate phase to 
a homogeneous density phase as p is increased beyond a critical value p c . A mean-field 
approximation does not predict a phase transition and describes only the condensate 
phase. We provide heuristic arguments for understanding the numerical results. 
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1. Introduction 

One dimensional lattice models are important tools in nonequilibrium physics. They 
are used to model physical systems like shaken granular gases [TJ [2], vehicular traf- 
fic [3] and ribosomal motion on mRNA [4j. The scope of these models however is larger 
than modelling particular phenomena. In absence of a general formalism, their study 
helps us understand the nature of steady states and phase transitions in nonequilib- 
rium systems [5]. In particular, the asymmetric simple exclusion process (ASEP) [6] 
and the zero range process (ZRP) [7] have emerged as paradigms and have provided 
us important insights owing to their simplicity and analytic tractability. The ASEP 
with open boundaries shows the interesting phenomenon of boundary induced phase 
transitions [8]. The phase behaviour and steady state properties in this case can be 
obtained analytically [9j. Apart from open boundaries, dynamical features also have a 
significant effect on the steady state. The steady state of the ZRP is known to be a 
product measure for certain forms of the hop rate, allowing analytic treatment j2j HP] . 
Other dynamical features like chipping and aggregation [HI [121 H31 El> disorder [15] 
and evaporation/deposition [16] have been shown to affect the steady state behaviour. 
Our work involves studying the effect of non-local hops, a feature whose effect on simple 
exclusion process has recently been explored [T7J HH [191 (20J [21]. The study of ZRP like 
models with hops of arbitrary length has been done in the equilibrium context in the 
urn models [22J. 

In this paper we study a one dimensional lattice model similar to the ZRP, allow- 
ing for unlimited occupancy. The particles have an attractive interaction with hopping 
rates inversely proportional to the occupancy rii, i being the label of the site to which 
they belong. The dynamics consists of two possible moves - a particle from a randomly 
chosen site can either move to one of the nearest sites with a probability (1 — p)/2 or 
make a non-local hop to one of the two sites at a distance of with probability p/2 (see 
Fig. (CQ)). Thus the non-local move is always to a distance equal to the occupancy of the 
chosen site. While the attractive interaction aids the clustering of particles on a site, 
the non-local hops promote declustering. We show that the competition between these 
opposing mechanisms leads to a phase transition between a condensate and a homoge- 
neous density phase. It is well known for the ZRP that if the hop rate decays to zero 
in the limit of infinite occupancy, there is a condensate formations at all densities [2]. 
Thus the p — > limit of our model has a steady state with a condensate containing a 
finite number of particles, while the density on the rest of the lattice approaches zero. 
As we will see, the steady state changes as p is increased and the condensate disappears 
at a finite value of p. 

Our model can be compared to the conserved mass aggregation (CA) models stud- 
ied in [HI [12l [13] where a phase transition occurs due to the opposing mechanisms of 
chipping and aggregation. In these models, a cluster of particles moves as a whole and 
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merges with particles already present at the arrival site. This promotes clustering. On 
the other hand, it is also possible for single particles to chip off from a cluster, promot- 
ing declustering. The control parameters here are the density and the probability of 
chipping. As the density is increased for a given chipping probability, the system goes 
from a homogeneous to a condensate phase. 

Our results show some similarity with the model described above. We find a phase 
transition controlled by two parameters - the density and the probability of non-local 
hops (Phase diagram, Fig. (j3j)). While nonequilibrium models with a condensation 
transition have been well studied, the importance of this work lies in a novel physical 
mechanism for the transition. The melting of the condensate in our case crucially de- 
pends on the symmetry in particle dynamics and on the dimension. The phase transition 
happens only with symmetric dynamics and is expected to occur only in one dimension. 
We will describe the physical mechanism behind these results in the discussion section. 

The paper is organised as follows. In section 2, we describe the model and present 
numerical results for the phase transition. We analyse these results in section 3. We 
find that a mean field approximation for the model does not capture the phase transi- 
tion. We provide heuristic arguments for the physical mechanism of the transition. Our 
arguments are supported by numerics. We summarise our results in section 4. 

2. Description of the model, Numerical results 

2.1. The model 

We consider a lattice with L sites and iV particles moving on it. We take periodic 
boundary conditions and allow for arbitrary occupancy at the sites. The dynamics is as 
follows: we choose a site at random and if it is occupied, we attempt to move a particle 
from this site. The rate of a particle moving out of this site, u(n), is determined by the 
occupancy n of the site and we take u(n) = 1/n. A particle that leaves the site can 
either make a local move to one of the nearest neighbour sites with probability (1— p)/2, 
or it can make a non-local move to one of the two sites at a distance n with probability 
p/2. Here p lies in the interval [0, 1]. The dynamical moves with the overall rates of 
particle motion are shown schematically in Fig. [TJ The non-local hops and periodic 
boundary condition impose a limit on the number of particles. Long jumps of length 
L/2 or more will move the particle closer to the originating site instead of taking it 
farther. Therefore, the maximum possible value of N is ^ — 1. 



Condensation transition in a model with attractive particles and non-local hops 
u(n)p/2^ ^ ^ ^u(n)p/2 



o 
o 
o 
o 
o 



o 



u(n)( 1 -p)/2/^ \u(n)( 1 -p)/2 

o o 
o o o 



o 
o 



Figure 1. Schematic diagram showing the possible dynamical moves. A particle 
(shown in grey) from a randomly chosen site can attempt one of the four possible moves 
shown with the corresponding probabilities. The value of n in the above example is 4. 



2.2. Numerical results 

We performed Monte-Carlo simulations using the dynamics described in the previous 
section. We computed the distribution P(n) which gives the probability of finding n 
particles on any given lattice site. Fig. (T5]) shows the variation of P(n) as a function of 
the parameter p. The number of particles here is N = -| — 1 (L = 2048, N = 1023) . For 
smaller values of p we see a peak in the distribution function at large n, in addition to a 
continuous function at small n. The peak represents a condensate while the continuous 
part represents the distribution of particles which are not part of the condensate. The 
peak disappears as p is increased beyond a critical value p c indicating that there is no 
condensate for p > p c . The value of p c in this case is 0.24. The probability distribution 
at p c is a power law and the exponent —2.3 ±0.1 seems to best fit the data at small ar- 
guments. The data for p > p c does not fit a simple functional form. The noncondensate 
part of the probability distribution for p < p c also does not fit a simple function unlike 
previously studied models where it is shown to be a power law [TU [12] . 



Figure ([3]) plots the average number of particles in the condensate, N c , as a function 
of the system size L, for density close to half (N = ~ — 1). The data fits a straight 
line indicating that the number of particles in the condensate is a finite fraction of the 
total number. This fraction can be inferred from the slope of the lines. As expected 
from the ZRP results, at p = nearly all of the mass is contained in the condensate. As 
p is increased, the fraction of particles in the condensate decreases but is still a finite 
fraction of the total mass. Thus the picture of the steady state in the < p < p c regime 
is that in addition to a condensate, there is a finite number of particles distributed 
homogeneously on the lattice. 



We computed the probability distribution P(n) for different values of particle 
density p (= N/L) and studied the change in steady state as a function p. Our results 
show that there is a phase transition and the value of p c decreases with the decrease in 
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Figure 2. Probability distribution P(n) versus n for (a) (top left) p < p c , with p = 0.0 
(red), 0.1 (green), 0.2 (blue), (b) (top right) p = p c where p c — 0.24 and the straight 
line shows a power law with exponent —2.3 and (c) (bottom) p > p c with p = 0.5 
(red), 1.0 (green). 
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Figure 3. Average number of particles in the condensate, N c , plotted versus the 
system size L for p = 0(red), O.l(green), 0.16 (blue). The total number of particles is 
N = -§ — 1. The slope of the straight lines indicates that the fraction of mass contained 
in the condensate, ^f-, is 1, 0.86 and 0.7 for p — 0, 0.1 and 0.16 respectively. 
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Figure 4. The figure shows the phase boundary between the condensate and the 
homogeneous density phases. The x-axis shows the density p = j- and the y — axis 
shows p. 

density. The phase diagram in Fig. (j4j) sums up our results. While we have kept the 
density constant and studied the phase change as a function of p, one could equivalently 
see the transition from a homogeneous density to a condensate phase by keeping p 
constant and increasing the density. 

3. Discussion 

The factorisation of steady state probability in the zero range process makes analytic 
treatment possible. The criterion for factorisability has been discussed in [2] and it has 
been shown that the steady state factorises if the hop rate from a site / to a site k, 
Uki(ni), has the general form Uki(ni) = Ui{ni)Wu- Here Ui(n{) gives the probability of a 
particle leaving the site I which is occupied by n\ particles and Wu is the probability 
that the particle hops to site k. In our case however, there is an explicit dependence of 
Wki on ni and the steady state in our model can not be written in a factorised form. 
When exact treatment is not possible, mean-field analysis has proven to be successful 
in studying phase transitions in similar systems [HI Ej. We will show that in our case, 
a mean-field equation that assumes no correlation between the sites does not capture 
the phase transition and the steady state behaviour. We will provide a heuristic under- 
standing of the phase transition mechanism and discuss a possible cause of the failure 
of the mean-field treatment. 

We can write the equation for the probability P(n) in the mean-field approximation 

as 
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dP(n) ( N ~ n N ~ n \ 

= - P{n)u{n) - P{n) (1 - p) £ P{k)u{k) + p £ P{k)u{k) (1) 

\ fc=i fc=i / 

/ N—n N-n \ 

+ P(n + l)«(n + 1) + P(n - 1) (1 - p) ^ P{k)u{k) + p J2 P{k)u(k) 

\ k=l k=l / 

Here = 1/n is the probability of a particle escaping a site having n particles. 
The probability P(n) will decrease if a site occupied with n particles either loses or 
gains a particle. The first term on the RHS is a loss term arising from a particle leaving 
a site having n particles. The second term accounts for a particle being added to the 
n particles already present. This may happen in two ways; via a local hop from the 
nearest neighbours or from a non-local hop from any site on the lattice. Thus the term 
with co-efficient (1 — p) inside the bracket comes from the local jumps and the sum- 
mation is over all possible occupancies for the nearest neighbour site. The term inside 
the bracket with co-efficient p arises from non-local hops, the summation here is over 
sites at a distance k which have exactly k particles. The third and fourth terms are 
gain terms - P(n) increases if one of the n + 1 particles present at a site leaves or if a 
particle enters a site having n — 1 particles. The last term again has local and non-local 
contributions. In all terms, there is an implied left and right contribution of 1/2 and 
1/2 each, which add up to give 1. 



The above equation reduces to 



dP(n) N ~ n 

—±-± = - P( n )u{n) - P{n) £ P{k)u{k) + P{n + l)u{n + 1) 
dt k =i 

N-n 

+ P{n - 1) J2 P{k)u{k) (2) 

k=l 

which is independent of p and carries no information about the non-local dynamics. 
In fact Eq. (|2J) is the same as that for a simple ZRP with nearest neighbour hops i.e. 
the p=0 limit of our system. It therefore describes a steady state with a condensate 
containing all the particles. Thus, the above equation does not describe phase transition 
in the one dimensional system we are looking at. Since mean field theory is expected 
to work in higher dimensions, the observed transition should be a feature restricted to 
low dimensions. 



A comparison of the relevant time scales in a system can provide insight into the 
cause of the phase transition [13] and we employ a similar method to analyse our system. 
Consider the system to have a condensate of mass n' at some value of p = p'. Since the 
rate of a particle to leave the condensate is u{n') = 1/n', one particle is lost every n' 
time steps. A fraction p of these particles make non-local hops, therefore there is one 
non-local hop in n'/p' time steps. Now consider the time of return to the condensate 
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of the particle that has made a non-local hop. Since the density on the lattice is small, 
this particle will return diffusively, making mainly local hops. Thus the time of return 
of this particle goes as n' 2 . This is comparable to the time it would take for the whole 
condensate to melt down, t' = ^rY^i'i ~ n ' 2 - The local hops are irrelevant as the 
particles can return to the condensate in a time of order one and therefore cannot cause 
declustering. 

So we see that it is an interplay of the two time scales of condensate melting versus 
diffusive return that determines the phase of the system. Which one of these domi- 
nates is determined by the value of p as well as the density of particles on the lattice. 
A higher value of p leads to more declustering and a smaller density leads to smaller 
condensates, both of which favour a homogeneous density phase. On the other hand a 
higher density and a smaller value of p will favour condensate formation. This heuristic 
argument, if true, leads to interesting consequences which serve as checks for its validity. 

Suppose the probability of a particle leaving a site has a general form u(n) = l/n a , 
with our model above corresponding to the special case of a = 1. The argument of 
competing time scales predicts that for a > 1, the steady state will have a condensate 
for any value of p. This is borne out in our simulations, Fig. (J5J). For a < 1, the p = 
state is a condensate phase as predicted by the ZRP but we expect that there should 
be no condensate for p > 0. We indeed see that for small a, the condensate disappears 
at very small values of p, Fig. (JHJ) (inset). Another possible alteration in our original 
model can be that the non-local hop, instead of having a length n, has a length equal 
to the integer part of n 7 . For 7 < 1, we expect from the argument of time scales that 
the condensate will always persist. We have also confirmed this prediction with our 
numerics. These results give credence to our argument. 

We also studied the effect of introducing asymmetry in the dynamics. In the model 
above, the particles have an equal probability of moving either in the left or the right 
direction. We introduced an asymmetry so that the probability of moving in one direc- 
tion is greater than the other. Because of this asymmetry, particles making a non-local 
hop can return to the cluster in a time proportional to size of the cluster n' instead of 
the diffusive n' 2 . Thus the asymmetry acts to stabilise the condensate. Our results in 
Fig. (EJ) show that the introduction of even a slight asymmetry causes the homogeneous 
density phase to disappear. We again see that it is the diffusive dynamics which is 
responsible for the phase transition. 

Let us now consider dimensions higher than one. Since we are working with a finite 
density, the number of particles on the lattice scales as L d , d being the dimension. Thus 
the size of the condensate also scales as L d . On the other hand, our dynamical rules 
do not allow the particles to make jumps of length larger than the system size L. Thus 
a particle farthest from the condensate will make a diffusive return in time of order 
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Figure 5. The main figure shows P{n) versus n for a = 1.5 and p — 0.0 (red), 0.5 
(green), 1.0 (blue). The condensate persists for all values of p. The inset shows P(n) 
versus n for a = 0.5 and p — 0.0 (red), 0.03 (green), 1.0 (blue). There is a condensate 
phase for p = but it vanishes at a very small value of p. 

L 2 . This is always smaller than the time scale of cluster melting which goes faster than 
L d+1 . Thus, if we define a restricted model in higher dimension, it will always have a 
condensate as predicted by the mean field theory. Our model therefore is expected to 
work only in one dimension. 

4. Conclusions 

We have studied phase transition in a nonequilibrium lattice model in one dimension. 
The transition from a condensate phase to a homogeneous density phase is caused by 
two opposing factors - the attractive interaction between the particles which promotes 
condensation and the non-local move which destabilises the condensate. We find that the 
diffusive dynamics is crucial to the transition, as is the fact that the leaving probability 
u(n) is inversely proportional to the number of particles in the cluster, u(n) = 1/n. 
This explains why the transition takes place only in one dimension and the mean-field 
treatment fails. For the same reason, an introduction of asymmetry in the dynamics 
leads to only one phase, the condensate phase. 
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Figure 6. P(n) versus n for a model with asymmetric dynamics. The particles 
preferably move to the right (probability 0.6) rather than symmetrically (probability 
0.5 on each side.). We show data for p = 0.0 (red), 0.5 (green) and 1.0 (blue). We 
see that the condensate exists for all values of p. The system size L = 1024 and the 
number of particles N = 511. 
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